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Gravitational Microlensing - A Report on the MACHO Project 

Will Sutherland 

Astrophysics, Dept, of Physics, University of Oxford, 

Oxford 0X1 3RH, U.K. 

There is abundant evidence that the mass of the Universe is dominated by dark matter of 
unknown form. The MACHO project is one of several teams searching for the dark matter around 
onr Galaxy in the form of Massive Compact Halo Objects (MACHOs). If a compact object passes 
very close to the line of sight to a background star, the gravitational deflection of light causes 
an apparent brightening of the star, i.e. a gravitational ‘microlensing’ event. Such events will 
be very rare, so millions of stars must be monitored for many years. We describe our search for 
microlensing using a very large CCD camera on the dedicated 1.27m telescope at Mt. Stromlo, 
Australia: currently some 14 events have been discovered towards the Large Magellanic Cloud. The 
lack of short-timescale events excludes planetary mass MACHOs as a major contributor to the dark 
matter, but the observed long events (durations 1-6 months) suggest that a major fraction may 
be in fairly massive objects ~ 0.5 M©. It is currently difficult but not impossible to explain these 
events by other lens populations; we discuss some prospects for clarifying the nature of the lenses. 
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I. INTRODUCTION: DARK MATTER 


This section gives an overview of the evidence for dark matter. This is a very large subject so only a brief outline 
can be given here. Further details and references may be found in, e.g., Peebles ( 1993 ). 

There are several strong lines of observational evidence for the existence of large quantities of dark matter in the 
universe. This is often parametrised in units of the critical density by O = p/pcrit, where p is the average density of 
matter in the universe, and pcrit = ^Hq/SttG = 2.78 x 10^^ Mq Mpc“^ is the critical density at which (for zero 
cosmological constant) the Universe is balanced between indefinite expansion and eventual recollapseQ. Generally Oq 
denotes the total matter density, and, e.g., flB^^sturs denote the fraction of critical density contributed by baryons, 
stars, etc. 


A. Evidence for Dark Matter 


The rotation velocities of spiral galaxies as a function of galactocentric distance can be accurately measured from 
the Doppler effect: at large radii where the stellar surface brightness is falling exponentially, velocities are obtained 
for clouds of neutral hydrogen using the 21 cm hyperfine line. The resulting ‘rotation curves’ are found to be roughly 
flat out to the maximum observed radii ^ 30kpc, which implies an enclosed mass increasing linearly with radius. 
This mass profile is much more extended than the distribution of starlight, which typically converges within ^ lOkpc; 
thus, the galaxies are presumed to be surrounded by extended “halos” of dark matter (e.g. Ashman 1992). 

Perhaps the most compelling evidence for dark matter comes from clusters of galaxies. These are structures of 
~ 1 Mpc size containing ^ 100 galaxies, representing an overdensity of > 1000 relative to the mean galaxy density. 
They may be assumed to be gravitationally bound since the crossing times for galaxies to cross the cluster are only 
^ 10% of the age of the Universe. Their masses can be estimated in three independent ways: 

i) From the virial theorem using the radial velocities of individual galaxies as ‘test particles’. 

ii) From observations of hot gas at ^ lO^iF contained in the clusters, which is observed in X-rays via thermal 

bremsstrahlung. The gas temperature is derived from the X-ray spectrum, and the density profile from the 
map of the X-ray surface brightness. Assuming the gas is pressure-supported against the gravitational potential 
leads to a mass profile for the cluster. 


iii) From gravitational leasing of background objects by the cluster potential. There are two regimes: the ‘strong 
lensing’ regime at small radii, which leads to arcs and multiple images, and the ‘weak leasing’ regime at large 
radii, which causes background galaxies to be preferentially stretched in the tangential direction. 


All these methods lead to roughly consistent estimates for cluster masses, (e.g. Carlberg et al. 1998, Blandford & 
Narayan 1992); visible stars contribute only a few percent of the observed mass, and the hot X-ray gas only ^ 10—20%, 
so the clusters must be dominated by dark matter. 

On the largest scales, there is further evidence for dark matter: ‘streaming motions’ of galaxies (e.g. towards nearby 
superclusters such as the “Great Attractor” ) can be compared to maps of the galaxy density from redshift surveys to 
yield estimates of ( $trauss fc Willick , 1995 ). Here the theory is more straightforward since the density perturbations 
are still in the linear regime, but the observations are less secure. A similar estimate may be derived by comparing 
our Galaxy’s 600kms“^ motion, measured from the temperature dipole in the cosmic microwave background (CMB), 
to the dipole in the density of galaxies. 

There are also some useful guidelines from theory. Primordial nucleosynthesis successfully explains the abundances 
of the light elements '^He, D, ^He & "^Li if the density of baryons satishes Hb ~ (0.01 — 0.05) (/i/0.7)“^. This suggests 
that baryons do not dominate the universe, but (depending on the controversial D abundance) this is probably higher 
than the density of visible matter Dvis ~ 0.01, so allows the dark matter in galactic halos to be mainly baryonic. 

Furthermore, it is easier to reconcile the observed large-scale structure in the galaxy distribution with the smallness 
of the microwave background anisotropies if the universe is dominated by non-baryonic dark matter. The theory of 


^We use standard astrophysical units throughout: G is Newton’s constant, M© « 2.0 x 10®° kg is the mass of the Sun, 1 AU 
~ 1.5 X 10® km is the mean Earth-Sun distance, 1 parsec (pc) = 648,000/7r AU, 3600 arcsec = 1 degree. Ho is the Hubble 
constant, parametrised hy h = Ho/{100 kms“® Mpc”*^) ~ 0.5 — 1.0. Note that the Sun is « 8 kpc from the center of the Galaxy, 
the Magellanic Clouds are at « 50 kpc, and the Andromeda galaxy at « 800 kpc. 
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inflation (postulated to solve the horizon and flatness problems) prefers a flat universe with ffg + = 1, where is 

the dimensionless cosmological constant; thus flo = 1 is the most ‘natural’ value, which seems to require non-baryonic 
dark matter. The predictions of inflation should be testable in the next decade with observations of CMB anisotropy 
by the MAP and Planck Surveyor satellites. 


B. Dark Matter Candidates 


Some ‘obvious’ dark matter candidates are excluded by a variety of arguments (Carr, 1994): hot gas is excluded by 
limits on the Compton distortion of the blackbody CMB spectrum; atomic hydrogen is excluded by 21 cm observations; 
ordinary stars are excluded by faint star counts; ‘rocks’ are very unlikely since stars do not process hydrogen into 
heavy elements very efficiently; hydrogen ‘snowballs’ should evaporate or lead to excessive cratering on the Moon; 
and black holes more massive than ^ 10® Mq would destroy small globular clusters by tidal effects. 

Most viable dark matter candidates fall into two broad classes: astrophysical size objects called Massive Compact 
Halo Objects or MACHOs, and subatomic particles (a subset of which are called Weakly Interacting Massive Particles 
or WIMPs). 

Each of these classes contains various sub-classes: for MACHOs, the most obvious possibility is substellar Jupiter- 
like objects of hydrogen and helium less massive than 0.08 Mq. Below this limit, the central temperature never becomes 
high enough to ignite hydrogen fusion, so the objects just radiate very weakly in the infrared due to gravitational 
contraction; thus they are usually known as ‘brown dwarfs’. Other MACHO candidates include stellar remnants such 
as old (and hence, cool) white dwarfs, neutron stars, and black holes (either primordial or remnants). 

For the particle candidates, they must clearly be weakly-interacting to have escaped detection, so possibilities 
include the ‘axion’ (hypothesised to solve the strong CP problem^), a neutrino (if one or more flavours has a mass 
~ 10 eV), and the popular ‘neutralino’ which is the lightest supersymmetric particle, thought to be stable. (Note 
that the term ‘WIMP’ is usual l y rese rved for the latter particle). There are active searches in progress for all of these 
particles (e.g. Jungman et al. ( 1996| )), but they will not be discussed here. 


II. GRAVITATIONAL MICROLENSING 


Even if MACHOs comprise most of the Galactic dark matter, they will be very hard to detect directly since 
they would emit very little electromagnetic radiation; future infrared searches may be able to constrain part of the 
parameter space, but not all. Thus, it is more promising to detect their gravitational field, via its influence on the 
light from background sources. In addition to the well-known test of General Relativity (GR) by light deflection by 
the Sun, gravitational leasing now has man y applications in cosmology: the first example of a quasar doubly imaged 
by an intervening galaxy was discovered by Walsh, Garswell fc Weymaim| ( 1979 ), and some 20 such objects are now 
known. More recently there have been many discoveries, many using the Hubble Space Telescope, of lensing of distant 
galaxies by intervening clusters. This takes various forms: sometimes multiple images are observed, in other cases 
highly distorted ‘giant arcs’ are found, while weaker image distortions are found at larger separations as discussed in 
Section ^ An overview of recent gravitational lensing observations is given by Kochanek & Hewitt (1996). 

The principle of lensing by MAGHOs in our Galaxy is very similar, but we shall see that the relevant angular 
separation is much smaller so the observable consequences are quite different. 


A. Principle of Microlensing 

If a compact object of mass M at distance I lies exactly on the line of sight to a (small) background source at 
distance L, the light deflection by GR causes the source to appear as an ‘Einstein ring’ (cf Hewitt et a71988) with 
‘Einstein radius’ te in the lens plane; the (small) light deflection angle is a = AGM/c^te, and geometrical optics 
gives a = teII -f rEl{L — 1), thus 


^ The ‘strong CP problem’ is that CP violation in the strong interaction is very small. Limits on the neutron electric dipole 
moment require an arbitrary QCD phase angle to be zero within 1 part in 10®. This seems unlikely by chance; it is arranged 
by the ‘Peccei-Quinn mechanism’ leading to the axion. 
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where x = l/L is the ratio of the lens and source distances; the corresponding Einstein angle is 9e = te/I- [Note that 
te \/rsL, where rs is the Schwarzchild radius of the lens.] As we introduce a small misalignment angle /3 between 
lens and source, it is clear that two images will be formed on opposite sides of the lens, collinear with the lens and 
source, at angular positions 

0± = O.5 (^/3±^/32 + 40|^ (2) 


from the lens. 

For a source star at SOkpc and a lens at lOkpc, the Einstein radius is xe ~ 8^/M/ Mq AU, and 6e ~ 10“^ 
arcsec. This is far below the resolution of ground-based telescopes, for which O^es ~ 1 arcsec (5p,rad) which is set by 
atmospheric turbulence or ‘seeing’. Thus the doubling of the star’s image is not observable, hence the general term 
‘microlensing’. However, leasing preserves surface brightness, thus the apparent flux of the source is magnified by the 
ratio of the (sum of the) image areas to the source area. A typical star at 50 kpc has an angular radius Os ^ 10“® 
arcsec, and thus may usually be treated as a point source. The fact that 

Ores 0e'> Os (3) 


explains much of the simplicity of microlensing; the first inequality gives the ‘micro’lensing, while the second inequality 
gives the simple magnification formula of Eq. ^ below. 

For a point source and any axisymmetric lens, the magnification A of each image is simply 


Oi dO 

/3i dp 


For the point lens, the total magnification factor is (Refsdal, 1964) 


A = A+ + 


vP + 2 

u'PvP + 4 


( 4 ) 

( 5 ) 


where u = P/Oe = O/te is the misalignment in units of the Einstein radius, and h is the distance of the lens from the 
undeflected line-of-sight^. 

Equation || gives A « for u ^ 0.5 and A k. 1 -\- 2u~^ for u ^ 2; thus the magnification may be very large, 
but is only appreciable for u ^2. Of course a constant magnification is not usually measurable, but since MACHOs 
must be in motion in the gravitational potential of the Galaxy, the magnification will be time-dependent due to the 
changing alignment; thus a microlensing ‘event’ will appear as a transient brightening with a timescale 

t = ~ 140\/M/ M0 days; (6) 

VE 

where v± ~ 200kms“^ is the transverse velocity of the lens^ relative to the (moving Earth-source line. Assuming 
constant velocities, the apparent ‘lightcurve’ A(t) of the event is simply given by Eq. @ with 

u{t) = [u^i„ -k [2{t - tmax)/t^^] (7) 


where the minimum misalignment Umin (thus maximum magnifica tion) occur s at t ime tmax- This application of 
microlensing to probe Galactic dark matter was first suggested by Paczynsk: ( 1986| ), so this form of A{t) is often 
known as the ‘Paczynski curve’. 

An example of the (unobservable) image behaviour for an event with Umin = 0.15 is shown in Fig. |^, and the 
lightcurves A(t) for various values of Umin are illustrated in Fig. ^ Note that for large magnifications, the events have 
a distinctive shape with a sharp central peak and broad wings. 


® Two points are noteworthy. Though A —> oo as w ^ 0, the average magnification integrated over a finite source is 
well-defined. Also, though A > 1 for all u, this does not violate energy conservation since introducing the lens changes the 
background metric. 

^This is estimated from the circular velocity of the Galaxy. Though the orbits of dark matter objects are probably not circular, 
their typical speeds must be of this order by the virial theorem. 
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FIG. 1. A microlensing event seen at ‘perfect’ resolution. Axes show angular offsets on the sky of the source from the lens 
(dot) in units of the Einstein angle 9e (defined following Eq. ^); the dashed circle is the Einstein ring. The series of open circles 
shows the ‘true’ source position at successive timesteps, for an event with Umin = 0.15. Eor each source position, there are two 
images (solid regions) collinear with the lens and source, as indicated by the straight line. 



Time (r,, / Vj^) 


FIG. 2. Lightcurves A{t) of microlensing events for six values of the impact parameter Umin = 0.0, 0.2,..., 1.0 as labelled. 
Time is in units of Einstein radius crossing time veIve = t /2. The inset illustrates the Einstein ring (dashed) and the source 
paths relative to the lens (dot) for the 6 curves. 


6 


















The timescale t of Eq. is observationally convenient, and the dependence oc '/M means that a large window 
in MACHO mass is accessible. This window is set at the low mass end by the finite size of the source stars, which 
limits the maximum magnihcation to « (1 + so lenses of mass ^ cannot produce appreciable 

magnifications. At the high mass end the limit ^ 100 Mq is set both by the patience of the observers and the falling 
event rate (see below). This range covers most of the plausible MACHO candidates. 

The microlensing ‘optical depth’ r for a given source is defined as the mean number of lenses within their own 
Einstein radius of the observer-source line; for r 1 as here, at most one lens gives an appreciable effect, so r is just 
the probability that a random star is microlensed with it < 1 or A > 1.34 at a given instant; thus 


T = / dl 


dM n{l, M) TTr‘^{l, M) 


AttGL 


dl 


dM M n{l, M) a:(l — x) 


4ttGL 


f dl p{l) a;(l — x) 
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( 8 ) 

(9) 

( 10 ) 


where n{l, M) dM is the number density of compact objects at distance I in the mass interval (M, M + dM), and p{l) 
is their mass density. Since xe oc VM, the ‘cross section’ of a lens at given I is oc M. Thus r depends only on the 
mass density profile p(l), and not on the individual lens masses. Using the virial theorem, it is strai ghtforward to 
show that T ~ ~ 10 “® where v is the orbital velocity of the Galaxy; more detailed calculations (Griest, 1991 ) 

using a realistic dark matter profile give 


tlmc « 5 X 10 


-7 


( 11 ) 


for microlensing of stars in the Large Magellanic Gloud by a ‘standard’ dark halo made entirely of MACHOs. 

Note that the rate of microlensing events T does depend on the lens masses via the durations; clearly the product 
of the event rate and the mean duration is proportional to the optical depth. We have 

r{t) = -r , ( 12 ) 

TT 

where the geometrical factor of 4/7r arises because t is defined as the time for the lens to move by one Einstein 
diameter, rather than the time it spends within the Einstein disk. For the LMC this leads to 

r « 1.6 X 10“®-\/ Mq/M eventsperstarperyear. (13) 


for a dark halo comprised entirely of MACHOs with mass M. Thus, high-mass MACHOs would produce very rare 
long-lasting events, while low-mass MACHOs would produce relatively more short-duration events. 


B. Microlensing Signatures 

The microlensing optical depth r ^ 10“® is very much smaller than the fraction of intrinsic variable stars ^ 10“^. 
However, microlensing has many distinct signatures which are very different from previously known types of variable 
star: 

i) Since the optical depth is so small, only one microlensing event should be seen in any given star. 

ii) Gravitational leasing is independent of wavelength, so the star should not change colour during the event. 

iii) The events should be symmetrical in shape and described by Eqs. || and ^ . These have 3 free parameters, but 

the event ‘shape’ only depends on Mmin, while t and tmax just represent a linear transformation of the time axis 
in Fig. |. 

In contrast, most variable stars are periodic or quasi-periodic, they are usually asymmetrical in time and show colour 
changes due to changing temperatures. 

Also, if many candidate microlensing events can be found, they should satisfy several statistical tests: 

iv) Microlensing does not discriminate between types of star, so the events should be distributed across the colour- 

magnitude diagram in proportion to the total number of stars. 
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v) The minimum impact parameter Umin should follow a uniform distribution between 0 and some experimental cutoff 

ut- This translates via Eq. ^ to a model-independent distribution in peak magnification Amax- 

vi) The peak magnification Amax and event duration t should be uncorrelated. 

In practice criteria (ii) and (iii) are slightly idealised and small deviations may be seen, while criteria (iv) - (vi) are 
modified by the experimental detection efficiency, but this can be accounted for as seen later. 

The well-specified shape of microlensing is useful for discriminating against variable stars, but has the drawback 
that it limits the information which can be extracted from each observed event. Of the 3 fit parameters Umin, t , tmaxi 
two give only the ‘uninteresting’ information of when and how close the lens approached the line of sight. All the 
desired unknowns M, l,v± of the lens are folded into the single observable t via Eqs. and |[ Thus we cannot uniquely 
determine the lens mass or distance from the lightcurve. If we assume a distribution function in distance and velocity 
for the lenses, the lens mass may be estimated statistically from the timescale, but with a large uncertainty; roughly 
a factor of 10^^ in M for a single event. Some methods for breaking this degeneracy are discussed later. 


III. OBSERVATIONS 

The very low optical depth above is the main difficulty of the experiment, and drives most of the observational 
requirements. 

Firstly, we require millions of stars at a distance large enough to give a good path through the dark halo, but 
small enough that the stars are not too faint. The most suitable targets are the Large and Small Magellanic Clouds 
(hereafter LMC and SMC), the largest of the Milky Way’s many satellite galaxies. Their location requires a Southern 
hemisphere observatory. The shape of the A(u) function means that it is optimal to monitor the maximum possible 
number of stars with ^ 10% photometric precision, rather than fewer stars with high accuracy. To get useful statistics, 
to monitor for long duration events and to check the uniqueness of candidates, frequent observations over several years 
are needed, thus a dedicated telescope is highly desirable; though a relatively modest (Im-class) telescope is sufficient. 
To check that the observed brightening is independent of wavelength, it is useful to have simultaneous observations 
in two different wavelength regions. 


A. The Telescope &: Cameras 


These requirements are incorporated in the MACHO experiment as follows: we have full-time use of the 1.27-metre 
telescope at Mt. Stromlo Observatory near Canberra, Australia from 1992-1999; it was brought out of mothballs and 
refurbished for this project. An optical corrector is installed near the prime focus to give good images over a wide held; 
this contains a dichroic beamsplitter to give simultane ous images in ‘blu e’ and ‘red’ passbands, covering wavelength 


ranges approximate ly 470 — 630nm and 630 — 760nm ( Hart et al. , 1996 ). Each focal plane is equipped with a very 
large CCD camera (Stubbs et al., 1993) containing 4 Loral CCD chips of 2048^ pixels; the held of view is a square of 


side 0.7 degrees. Each chip has 2 amplihers, so the images are read out through a 16-channel system, giving a readout 
time of 70 sec for 77 MB of data; the readout noise is 10e“ which is small compared to photon noise from the night 
sky. The typical exposure times are 300 sec for the LMC, 600 sec for the SMC and 150 sec for the Galactic bulge, so 
around 60-100 images are obtained per clear night; the stellar detection limit is around 21st visual magnitude, about 
10“® times fainter than the human eye. All the raw data are archived to Exabyte tape. The LMC subtends roughly 
6 degrees, so about 80 images are required to cover it. The observing strategy has varied so that some helds were 
monitored several times per night to be sensitive to short-duration events, while all helds are observed at least ~ 20 
times per year to maximise sensitivity to long-duration events. 


B. Photometry 


A typical night’s observing produces some 80 images containing up to 600,000 stars each, so high-speed software 
is required to measure their bri ghtnesses. We use a spec i al-pu rpose code called SoDOPHOT, customised from the 
well-known DOPHOT package (|Schechter, Mateo fc Saha , 1993 ). The major modihcation is the use of ‘templates’ as 
follows: for each held we choose one image from very good sky conditions as a ‘template image’. A ‘full’ reduction is 
run on this image to produce a list of detected stars, a ‘template’. Subsequent images are divided into 64 sub-frames 
called ‘chunks’ in each colour, to minimise point-spread function (PSF) variation and focal-plane distortions. For each 
chunk, ^ 30 bright reference stars are located and used to dehne a PSF, and a coordinate transformation and hux 
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scale relative to the template. Then, the brightness of all other stars is measured at their ‘known’ positions using the 
measured PSF, with neighbouring stars subtracted from the image. Each data point has 6 associated ‘quality flags’ 
such as the of the PSF fit and the fraction of the star’s image masked due to bad pixels, cosmic rays etc; these are 
used later to reject suspect data points. 

This use of ‘templates’ greatly speeds up the reductions since for each star there are only 2 free parameters, the 
flux and night sky brightness; a typical image reduction takes around 1 hour on a single CPU. We now process most 
of the image data to photometry measurements within 12 hours from the time of observation. This data is fed to a 
customised database on ~ 500 GB of RAID disks. 


C. Analysis 


The photometric time-series data is searched for variable stars and microlensing events; over 50,000 variable stars 
have been found, most of which are new. Thi s very large sample has many benefits for stellar physics, but is outside 


the scope of this article; see e.g. |Cook et al. (|l9971 ) for an overview. 


For the microlensing search, the lightcurves are first convolved with filters of various durations, and those showing a 
peak above some threshold are called ‘level-1’ candidates. These undergo a full 5-parameter fit to microlensing (where 
the free parameters are the 3 microlensing parameters Umin, t , tmax and the star’s red & blue baseline fluxes //{, /s), 
and numerous statistics are computed, including ‘total significance’ Ay^, goodness of fit in various windows, number 
of significant high points, etc. Stars satisfying weak cuts on these are output as ‘level-1.5’ candidates which can be 
browsed by eye, and then more stringent cuts are defined to select a set of final ‘level-2’ microlensing candidates. 
The definition of these cuts is necessarily somewhat subjective, since little was known in advance about the classes of 
variable star which bear some resemblance to microlensing; however, in practice it turns out that after eliminating two 
sparsely-populated regions of the colour-magnitude diagram (the upper main sequence and the very reddest stars), a 
uniqueness criterion requiring a constant baseline and a single high-significance brightnening with Amax ^1-5 appears 
to select a fairly clean set of microlensing events. 


D. Short History of Microlensing 


While this article focuses on results from the MACHO project, it is useful to revie w the progress of microlensing 
including a mention of the other projects; more details are given in Paczynski (1996). At the time of the proposal 
of Paczynski ( 1986| ), the project was not technically feasible, but following rapid advances in computers and CCD 
detectors, a seminar by C. Alcock at the Center for Particle Astrophysics, Berkeley in 1990 led to the formation of 
several teams. The first teams were MACHO, EROS and OGLE, and all three announced their first candidate events 


one towards the bulge by OGLE (Udalski et al. 
is much higher than towards the 


in late 1993: one event towards the LMC by MACHO (Alcock et al., 1993), two by EROS (Aubourg et al., 1993), and 

1993| ). It soon became clear that the event rate towards the bulge 
the world total is now around 150 bulge events and 15 LMC events. The 


.MC; 

majority of these events have been found by MACHO due to its dedicated telescope and larger data volume, and are 
discussed below; the OGLE and EROS teams have recently completed new dedicated t elescopes which wi ll increase 


the event rate. Several new groups have entered the field; DUO observing the bulge (Alard et al., 1996), AGAPE 


and VATT-Golumbia (Lomaney & Crotts, 1996|) observing the Andromeda galaxy, and MO A observing the LMC. 


Information on these projects may be found via the MACHO WWW 


IV. GALACTIC BULGE RESULTS 

It is a lucky coincidence that when the Magellanic Clouds are too low in the sky to observe, the dense ‘bulge’ of 
stars around the Galactic center is well placed; we spend about 1/3 of the observing time on the bulge. The line of 
sight to the bulge passes through the disk of our galaxy, where the mass density in stars is ~ lOx larger than the 
halo dark matter density. Detection of microlensing does not require that the lens be dark, just that it is not much 
brighter than the source star; so there must be a significant rate of leasing towards the bulge from ‘known’ low-mass 
stars (actually, these are not directly observed but must be present assuming a stellar mass function similar to the 


The MACHO WWW page is located at ittp://wwwmacho.mcmaster.ca 
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Time (days after 1993 January 12.5) 


FIG. 3. Lightcurves (observed flux vs time) for four representative microlensing events from the MACHO 1993 season bulge 
data. The data points with ±1 ct error bars show flux in linear units, normalised to the fitted baseline. For each event, the 
upper/ lower panels s how th e blue/red passbands. The smooth curve shows the microlensing fit, simultaneous to both colours. 
(From Alcock et al. (1997b)). 


solar neighbourhood). Thus, the bulge results are only indirectly relevant to the dark matter question, but are (a) a 
useful verification that microlensing really occurs and that our experiment can detect it, and (b) a probe of galactic 
structure and the low-mass end of the stellar mass function. 

At the time of writing, around 150 candidate microlensing events have been found towards the b ulge; 43 of these 


are from the first year’s data, for which a full statistical analysis has been done ( Alcock et ai , 1997b ). Some examples 
of these events are shown in Fig. and the Umin distribution shows excellent consistency with the predicted uniform 
distribution, after correction for the detection efficiency which is higher for small Wmin (large magnification). 


Over 100 additional events have been discovered by our ‘Alert System’ (Alcock et al, 1996a), whereby the photom¬ 


etry is carried out in real time within 12 hours of observation. Any star which is not in a pre-defined list of variables, 
shows a 7a upward excursion in both colours, and satisfies certain quality cuts is reported, and the full time series 
is extracted for human inspection. Events judged promising by eye are then announced as ‘alerts’, sent to an Email 


distribution list and placed on our WWW alert page at !ittp://darkstar.astro.washington.edu. 


The main conclusions from the bulge results are as follows: 

i) Microlensing is conclusively detected, since many of the events are of high quality and the statistical tests are well 
satisfied. Many events have been observed by multiple sites, and several have been observed spectroscopically 
in real-time and show no change in the spectrum (Benetti, Pasquini & West, 1995|) as expected. 


ii) Estimates of the optical depth towards the bulge range from 2.5 — 3.9 x 10 ° (Udalski et al, 1994a; Alcock ei 


at, 1997b) which is over twice as high as the previous theoretical predictions; this problem is alleviated if the 


bulge is actually a prolate ‘bar’ seen close to end-on, in agreement with recent dynamical evidence. 

iii) The distribution of event durations is roughly consistent with expectation assuming most of t he lenses are low-mas s 
stars. There is some evidence for an excess of short timescale events relative to predictions (Han fc Gould, 1996), 
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which may indicate a population of brown dwarfs in the bulge or may be due to blending (Sec. ^). 

iv) Some interesting events are found to show deviations from the ‘standard’ Paczynski curve; these are discussed in 
detail in Sec. M. 


V. MICROLENSING FINE STRUCTURE 

The standard ‘Paczynski curve’ assumes a single point source, single point lens and uniform relative motion between 
the lens and the observer-source line. These are good approximations for most events, but roughly 10 — 20% of 
observed events show significant deviations. At first glance one might think that these ‘deviant’ events cast doubt 
on the microlensing interpretation, but in practice many deviations occur for specific subsets of events and have 
been predicted theoretically, thus they actually help to prove microlensing. Deviations from the standard shape are 
especially valuable since they can provide additional information which breaks the intrinsic degeneracy between M, I 
and v± noted in Section ||. 


A. Blending 

The commonest but least interesting deviation is that caused by ‘blending’. Since the fields monitored are chosen 
to have a very high density of stars, there is a non-negligible probability that any observed ‘star’ actually consists of 
two or more stars within the observational seeing disk. Since 9e ^ 9res, usually only one of these will be lensed, so 
the observed magnification Aobs of the blend will be given by 

Aobs 

Aobs 1 

where /l, fu are the fluxes of lensed and unlensed stars within the seeing disk. If the colours of the lensed & unlensed 
stars are different, the event may appear chromatic, but there is still a linear relation between passbands. 

Blending can be estimated by extra fit parameters, but unfortunately there is a near degeneracy in that the lightcurve 
of a blended event appears very similar to that of an unblended event of lower A^ax and shorter t ; post-event Space 
Telescope images can help to break this degeneracy. Blending also affects the detection efficiency as discussed later. 


_ fu T Atrue/i 

fu + fu 

— {Atrue f) 


fu 


fu + fu 


(14) 


B. Parallax 


The Paczynski curve assumes uniform relative motion between the lens and the observer-source line. This is clearly 
inaccurate at some level due to the influence of the Earth’s orbit. For events in the observed region of parameter 
space (durations of months), the resulting deviation is small since a uniform component of the Earth’s motion is not 
separable from a change in velocity of the lens; only the change in the Earth’s velocity Auq over the duration of the 
event gives a measurable effect. The size of the effect scales approximately as Av(^/v, where v is the relative velocity 
of the lens ‘projected’ back to the solar system plane , v = u_l/(1 — x). 

This ‘parallax’ effect was predicted by Gould (|1992 ), and first observed by Alcock et al. ( 1995b ); this was the longest 
event in the MACHO year-1 bulge data, so such a deviation is not unexpected. The light curve of this event is shown 
in Figure The parallax effect is very useful since it enables the two components of v to be measured, by reference 
to the known parameters of the Earth’s orbit. Along with the usual t , one then has two constraints on the three 
unknowns M,l,v± of the lens; thus one obtains uniq ue relations for M, v± as a function of assumed lens distance 1. 
Additional constraints come from likelihood analysis ( Alcock et al , 1995b| ) since the probability of observing a given 
u is a sensitive function of 1. 


C. Finite Source Effects 


The usual approximation for A(u) assumes a point source. The full expression for a uniform circular disk is given 
by Witt & Mao (1994): in practice this is very close to the point-source formula except when f3 ^ 29s, be., the lens 
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FIG. 4. The lightcurve of a ‘parallax’ event; data shown as ±1 ct error bars. The dashed curve sh ows the ‘standard’ m l- 
crolensing fit of Eq. and the solid curve shows the fit including the effect of the Earth’s orbit. (From Alcock et al. (1995b)). 



Time (days after 1992 January 02) 


FIG. 5. The lightcurve of MAGHO Alert 95-BLG-30, which shows finite-source effects. Data points show flux relative to 
the fitted baseline with dblcr error bars. The upper p a rt of t he figure illustrates the geometry, with the source (small circle) 
and Einstein ring (dashed circle). (From Alcock et al. (1998c)). 
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comes close to transiting the disk of the source. Since Oe ^ for solar-mass lenses, this implies that finite-source 
effects should only be signihcant for very low-mass lenses or near the peak of high-magnification events. 


1998a); this was 


One very clear example of this effect has been observed, in MACHO Alert 95-30 ( |Alcock et al. 
an event with Amax ^ 25 involving a giant source star. The event was detected well before peak, and subsequent hts 
predicted a high peak magnihcation. Thus, the possibility of finite-source effects was anticipated in advance, and the 
peak was very frequently observed by some 5 observatories around the globe. 

The lightcurve data near the peak is shown in Figure P, which clearly shows the ‘shoulders’ due to the finite-source 
effect. This can be fitted with one additional free parameter, the angular size of the star in units of the Einstein angle, 
M* = 0s/0£;, which is 0.075 in this example. The apparent brightness and spectral type of the star give an estimate 
of Os', thus we obtain 9e, which leads to a mass-distance relation for the lens. 

Such finite source effects can be of considerable astrophysical interest due to the ‘differential magnihcation’ effect 
across the face of the star; thus, observations as the lens transits the source can constrain the center-to-limb variations 
in the star’s spectrum, which for giant stars is very hard to measure in any other way. 


D. Binary Lenses 


Undoubtedly the most spectacu l ar dev iation from the Paczynski curve is that arising from a binary lens (Schneider 
fc Weiss , 1986 ; Mao fc Paczynski , 1991 ). Although ^ 50% of all stars reside in binary systems, the binary leasing 
is most dramatic when the projected separation a is roughly comparable to the Einstein radius; so the expected 
frequency of observable binary effects for stellar lenses is reduced to ^ 10%. The hrst binar y lens discovered wa s 
OGLE-7 (Udalski et al., 1994b| ), which was soon confirmed independently in the MACHO data (Bennett et al, 1995). 

The binary lens is qualitatively different from a single lens since it contains astigmatism, which breaks the point 
singularity into one or more line singularities, i.e., ‘caustics’. These are closed curves in the source plane where the 
number of images changes by ±2. Caustics in the source plane map to ‘critical curves’ in the image plane where the 
determinant \d(3/d6\ = 0; hence the magnihcation is inhnite for a point source on a caustic. An example magnihcation 
map for a binary lens is shown in Figure the caustic is the 6-pointed closed curve. For a point source outside the 
caustic, there are 3 images; as the source crosses the caustic inwards, two new images appear ‘instantaneously’ with 
formally inhnite magnihcation, and these fade as the source is well inside the caustic. 



FIG. 6. Magnification map as a function of source position for a typical binary lens. Axes show sonrce offset from the lens 
centroid in units of Einstein angle Oe- Grey-scale shows magnification on a log scale, from 1 (black) to 10 (white). The lens 
consists of two point masses with 2 : 1 mass ratio separated by 19e (defined using the total mass of the system); the larger 
mass is at (-1/3, 0) and the smaller at (-1-2/3, 0). 

Clearly a binary lens can lead to a great diversity of lightcurves, since there are 3 additional free parameters; the 
mass ratio, the projected separation in units of the Einstein radius, and the orientation of the binary relative to the 
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source path. However, caustic crossings are common features, which must occur in inward-outward pairs, so these 
events are quite easily recog nised. Fitting of models to observed binary events is non-trivial due to the large parameter 
space and multiple minima (Mao & Di Stefanc, 1995). 


To date, about 8 definite binary-lens events have been observed, OGLE-7 above, DUO-2 (|Alard et alj [l996|) , 
MACHO LMC-9 below and about 5 by the MACHO alert system; there are several more probable cases where the 
data coverage is not quite sufficient to be certain. 

A related effect is the case of leasing by a star with a planetary system; here the caustics are much smaller, so such 
an event appears like a single-lens event with some probability (~ 1 — 20%) of a s hort-lived deviat i on if o ne image 
passes near a planet. This is valuable since it is sensitive to smaller planet masses (Bennett & Rhk, 1996) than the 
recent planet discoveries via radial velocity measurements. Thus, to detect planets via microlensing, it is desirable to 
take existing microlensing alerts and observe them much more frequently (~ hourly); this is being undertaken by two 
teams, PLANET and GMAN. 


VI. LMC RESULTS 


The microlensing searches towards the LMC and SMC are the most important for the dark matter question, since 
these lines of sight pass mainly through the outer Galaxy where the density is dominated by dark matter; thus, unlike 
the bulge case, the lensing rate for an all-MACHO halo is much larger than that from known stars. 

At present, we have analysed the first 2.3 years of data p] for our ‘high pri ority’ LMC fields; this sample contains 


some 8.5 million stars with 300-800 observations each (Alcock et al, 1997a). A set of automated selection criteria 


applied to this sample yields 12 objects; 2 of these are redundant detections of two stars which appear in field overlaps 
and are independently analysed, and a further 2 are rejected due to ‘magnification bias’ in that they were brighter 
than normal in our template image and have subsequently faded below our detection limit (one of these was almost 
certainly a supernova in a background galaxy). This leaves 8 candidate microlensing events, which are shown in Fig. |^. 
(These are numbered 1,4-10 since two low-quality candidates 2 & 3 appeared in an earlier paper but no longer pass 

the improved selection criteria). _ _ 

( 1993| ); the star’s spectrum shows no abnormalities 


Here event 1 is the event first announced in Alcock et al. 


(Della Valle, 1994). Event 4 was discovered by our Alert system, thus has more accurate follow-up observations 
(not shown here) which support the microlensing interpretation. Event 5 occurs in a faint star but has a very high 
peak magnification; the magnification appears slightly greater in the blue passband than in the red, probably due to 
blending. Event 9 is due to a binary lens; the fit shown in Fig. ^ is the single-lens fit used in event detection, thus is 
not appropriate. Events 6, 7 and 8 are good-quality candidates, while event 10 appears somewhat asymmetrical and 
may be a variable star; the inclusion or exclusion of this event has little effect on the results. The distribution of these 
candidates across the colour-magnitude diagram and in A^ax is consistent with the microlensing predictions given in 
Section |B. Some four of these candidates (1,4,5,9) are of high quality and are almost certainly due to microlensing; 
this suggests that most of the lower quality candidates are also microlensing, since if only the high-magnification 
candidates were actually microlensing while the others were variable stars, the implied distribution of Amax would be 
somewhat improbable. 

Given observed microlensing events, it is natural to ask whether the individual lenses can be detected post-hoc, 
either directly or by their microlensing other nearby stars. Unfortunately, direct detection is unlikely since the lenses 
may be 1 — 6 orders of magnitude fainter than the source star, and the relative angular velocity (“proper motion” in 
astronomical jargon) would be ^ 0.005 arcsec/yr; thus even after a decade they would be only marginally resolved 
by HST. Likewise, it will take centuries for a given lens to reach a neighbouring source star, and even then it is 
improbable that the alignment would be close enough to lens the second star. 

To obtain quantitative conclusions, we clearly need to assess our event detection efficiency. The detection probability 
for individual events is a complicated function of Amax, t , the observing strategy and the brightness of the source star, 
but all these distributions except t are known, and thus can be averaged over by a Monte-Carlo simulation. There 
are two levels of detail here: firstly a ‘sampling efficiency’ which assumes all stars are single sources and adds events 
at the time-series level; and secondly a ‘blending efficiency’ which adds artificial stars into the raw image data. The 
latter is more realistic, but typically only differs by ^ 15% from the sampling efficiency, since there is an approximate 


® Analysis of the 4-year LMC dataset is nearing completion, and yields about 6 additional events; there are also several further 
LMC events from the Alert system. These analyses and related efficiencies are not yet hnalised, but preliminary results indicate 
an optical depth similar or marginally lower than the 2-year value. 
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MACHO LMC 2-year Events 



0 200 400 600 

Time (days after 1993 January 12.5) 


FIG. 7. Lightcurves (observed flux vs time) for the eight candidate microlensing events from the 2-year Large Magellanic 
Cloud data. Data points (ilcr errors) show flux in linear units, relative to the fit baseline. For clarity, flux data have been 
averaged in time bins of width 1-4 days roughly scaled with the event durations (see labels) . Only data in one colour is shown. 
The curves show the single-lens microlensing fit to each event. (From Alcock et al. (1997a)). 


cancellation between the increased number of stars in a blend and the reduction in observed magnification. The event 
detection efficiency £{t) is shown in Fig. ||, and shows a broad peak between t ~ 30 — 300 days. 

Given a model for the halo density profile and velocity dispersion, it is straightforward to predict the event rate 
and distribution of timescales t for arbitrary lens mass; combined with the efficiency function, this gives an expected 
number of events iVexp- In practice it is convenient to show the expected number of events N{M) for an all-MACHO 
halo in which all MACHOs have a unique mass M, which is shown in Fig. Note that N{M) peaks at iV > 40 for 
M ^ 0.001 — 0.01 Mq; the reason is that for M ^ 0.1 Mq, most events are longer than 30 days where the efficiency 
curve is quite flat, while the event rate is falling oc For small masses M ^ 10“^ M©, the theoretical event rate 

is large but most events are shorter than 3 days where the efficiency is very low. The product of these two factors 
gives the peak in N{M). 

There are several complementary ways to analyse a set of detected microlensing events: firstly, one may obtain 
robust upper limits for given M by simply excluding models which predict too many detected events, shown in the 
lower panel of Figure 9. For substellar MACHOs we may obtain stronger limits since most events should have short 
duration {t < 20 days), but we have no candidate event shorter than this. The dotted lines in Fig. ^ show the expected 
number of events with t < 20 days, and the resulting limit: objects in the mass range 10“^ — 0.01 Mq comprise < 20% 
of a standard dark halo, and this applies whatever the shape of the mass function. We have extended this constraint 
to smaller masses using a separate search for very short duration events, called a ‘spike’ analysis; since many fields 
were observed twice per night, we can search for sets of four data points (two red-blue pairs) from a single night 


15 






































FIG. 8. Detection efficiency vs event duration S(t) for the LMC 2-year data. The dotted line shows the “photometric” 
efficiency assuming all stars are single, while the solid line incorporates the effects of blending. The efficiency is defined relative 
to the theoretical event rat e with itmin < 1, and thus contains a factor of 0.66 from the actual cut Amax > 1.75 or Mmin < 0.66. 
(From Alcock et al. (1997a)). 



10 -« 10-5 0.0001 0.001 0.01 0.1 1 10 


Macho mass (M^) 


FIG. 9. Upper panel: the expected number of microlensing events vs MAGHO mass N{M), including detection efficiency, 
for an all-MACHO dark halo with unique MACHO mass. The solid line shows total number of events, the dotted line shows 
events with duration t < 20 days. Lower panel: resulting limits on the fraction of the halo contributed by MACHOs of given 
mass; regions above the curves are exclud ed at 95% CL. The s olid line is derived from 8 events total, the dashed line from no 
observed events with t < 20 days. (From Alcock et al (1997a)). 
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which show a significant brightening in all four points. There is insufficient data here to obtain a microlensing fit, so 
if events were found we could not claim a detection; but using suitable selection criteria, no s uch events are found , 
which allows us to extend the excluded region down to around 10 ~^Mfn, as shown in F ig. [l0| (Alcock et ai, 1996b). 


Similar (nearly independent) limits have been obtained by EROS (Renault et al, 1997). 



FIG. 10. Upper limits on low-mass MACHOs based on the non-detection of very short duration ‘spike’ events, and no 
events with t < 20 days in the standard search. The lines show the 95% c.l. upper limit on the total mass of MACHOs within 
50 kpc of the galactic center, as a function of the assumed MACHO mass, for 8 different models of the dark halo. (From Alcock 
et o|] ( 1996b )). 


Secondly, one may obtain an unbiased estimate of the optical depth by 

ti 


"^est — 


AE 


E 


£{t^) 


(15) 


where ti is the timescale of the fth event and E is the ‘exposure’ in star-years. For the LMC 2-year sample this 


gives tlmc ~ 2.9lj!);g x 10“^. This quantity has the virtue of being independent of assumptions about lens masses 
and velocity distributions; the drawback is that it is subject to non-Poisson statistics, and also does not account for 
contributions from timescales outside the window of sensitivity where £ —^ 0. Thus, we caution that it is not valid 
to convert Test into a constraint on the total abundance of MACHOs without specifying an associated mass interval; 
unfortunately this over-general statement sometimes appears in the literature. 

The most dramatic result here is that the Test above is not much smaller than the value in Eqn. |^, and both 
Test and the number of 8 observed events substantially exceeds expectations from lensing by ‘known’ low-mass stars. 
Tensing by stars is estimated to contribute Tstars ~ 0.5 x 10“^ or around 1.1 event to the above sample; hence there 
appears to be a very significant excess. Thus, one may treat the events as a detection of halo dark matter and use 
a maximum likelihood model to estimate the most probable MACHO mass M and MACHO fraction of the halo /. 
Assuming a unique mass for all MACHOs, the result is shown in Fig. 
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and is M = 0.5^noMf7i 


f _ q c-t-O.S A 

-0.2 ^'■^0) / “ n.O_g 2. A 

is not simply the complement of 


common source of confusion here is that that the 95% c.l. excluded region in Fig. 

the 95% allowed region in Fig. |l^ The reason is that the unknown MACHO fraction of the halo is really a function 
iIj{M) (with f = J '4’iM)dM), not a point in 2-parameter space f,M. The likelihood analysis depends on the shape 
of ip{M) through the event timescales; this is unknown a priori and we must adopt some simple parametrisation, 
e.g. a (5-function or a truncated power law. If we assume a different shape for the mass function, the ‘allowed region’ 
in the likelihood analysis may change, but the excluded region in Fig. ^ does not. Thus one needs both analyses to 
extract the full information from the data. 
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FIG. 11. Likelihood contours for MACHO mass M and the fraction of the halo comprised of MACHOs /, for a standard 
model halo with a <5—function MACHO m ass distributi o n. Th e lower panel assumes all 8 events are due to halo lenses, the 
upper panel assumes 6 such events. (From Alcock et al. (1997E)). 


VII. DISCUSSION 


The interpretation of the LMC results is currently unclear. One conclusion is robust, that MACHOs in the planetary 
mass range 10“® to 0.05 Mq do not contribute a substantial fraction of the Galactic dark halo. Regarding the detected 
events, although the Poisson uncertainties are substantial, this is not the dominant uncertainty: the critical question 
is, to what population of objects do most of the lenses belong ? As noted above, the 8 events in the 2-year data are 
substantially in excess of the predicted ‘background’ of « 1.1 event arising from known stellar populations, which 
suggests that MACHOs in a dark halo are a natural explanation. 

However, there are some astrophysical difficulties with this interpretation, mainly arising from the estimated mass 
^ 0.5 Mq for the lenses. These cannot be hydrogen-burning sta r s in t he halo since such objects are limited to < 3% 
of the halo mass by deep star counts ( Gould, Bahcall fc Flynn , 1997 ). Modifying the halo model to slow down the 
lens velocities can reduce the implied lens mass somewhat, but probably not below the substellar limit 0.08 Mq. Old 
white dwarfs have about the right mass and can evade the direct-detection constraints, but it is difficult to form them 
with high efficiency, and there may be problems with overproduction of metals a nd overproduct i on of light at high 
redshifts from the luminous stars which were the progenitors of the white dwarfs (Chariot & Silk, 1995| ). Primordial 
black holes are a viable possibility, though one has to appeal to a coincidence to have them in a stellar mass range. 

Due to these difficulties of getting MACHOs in the inferred mass range without violating other constraints, there 
have been a number of suggestions for explaining the LMC events without recourse to a dark population: most of 
these suggestions construct some non-standard distribution of ‘ordinary’ stars along the LMC line of sight, to increase 
the stellar leasing rate. Fairly general arguments from faint star counts can limit the contribution from our Milky Way 

1997[). The contrib ution f rom L MC ‘self-lensing’, i.e. stars in the front of the LMC 


disk (Gould, Bahcall & Flynn 


leasing those in the back, is l ess wel l constr ained: Sahu (|1994| ) estimated an optical depth ^ 0.5 x 10“' from this, 
but this has been disputed by Gould ( 1995b| ) who obtains an upper limit of 0.25 x 10“^ using dynamical constraints. 
This can be tested given more data via the distribution of events on the sky, since lenses in the LMC itself should 
be concentrated towards the c enter o f the L MC. Other suggested lens populations includ e an unknow n dwarf galaxy 
roughly half-way to the LMC ( ^hao| , 1998a), a tidal tail of stars stripped from the LMC ( ^hac , 1998b), or a strongly 
flared and warped disk of our Galaxy (Evans et al, 1998). All of these proposals appear somewhat contrived but can 
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be tested observationally in the near future. 

However, the most decisive test for the lens population is to break the degeneracy in M, /, u_l and thus estimate the 
distance to at least a subsample of the lenses. This can be done for ‘non-standard’ lensing events using the methods 
outlined in Sec. but the percentage of events showing such deviations is expected to be small, especially for halo 
lenses. The ideal method for measuring lens locations is to discover lensing events in rea l time, as alrea dy happens, 
and then make extra observations from a small satellite at ^ 0.3 AU from the Earth ( |Gould , 1994b ). Since this 
distance is comparable to the Einstein radius of the lens, the event lightcurve will appear substantially different from 
the satellite and the Earth, e.g. the times of peak brightness are expected to differ by ~ 1 week. As with the parallax 
effect of Section this yields a measurement of the ‘projected’ velocity of the lens across the Solar system u_l/(1 — x). 
This can unambiguously determine whether the lens belongs to our disk, the halo or the LMC. Another possibility 
is that planned space astrometric missions su ch as SIM (P_aczynski , 199^ ) and GAIA should be able to measure the 
deflection of the light centroid ^ 10“"^ arcsec (Walker, 1995) during a microlensing event, which gives a measurement 
of the angular Einstein radius 9e and thus the lens angular velocity v±/l. Measuring both parallax and 9e for the 
same event gives a complete solution for M, I and v±. 

The prospects for microlensing appear bright: the MACHO project will continue to observe until at least 1999, 
and the EROS-2 and OGLE-2 telescopes have recently come on line, which should significantly increase the event 
rate, reducing the statistical uncertainties and enabling useful tests for the distribution of lenses across the LMC. 


(Alcock et al ., 

1997c; 

Palanque-Delabrouille et al.. 

1998) and candidate events towards M31 ( 

Tomaney & Crotts 

1996) 


The suggestions for stellar lensing populations can be tested observationally, though if these populations are not found 
we may need to wait a few years for one of the space-based missions to finally decide whether the lenses are in the 
dark halo. 
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